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ABSTRACT 

HARMS,  DEWEY  ELVIH.  Application  of  an  Objective  Analysis 
Scheme  to  Mesoscale  Observational  Network  Design.  (Under 
the  direction  of  DAVID  BARBER  and  GERALD  WATSON.) 

This  research  program  was  initiated  with  the  overall 
objective  of  estimating  the  probable  effectiveness  of  th^ 
proposed  surface  observational  network  C(PAM  , II )  which  will 
be  deployed  as  part  of  the  winter  1986  GALE  field  project. 
This  network  is  to  contribute  to  the  study  of  mesoscale 

features  within  developing  east  coast  cyclones,  and 

1 

especially  the  well-known  coastaL  front. 

To  accomplish  this  objective,  the  process  of  spatial 
weighted-averaging  interpolation  of  surface  data  was 
investigated  by  using  the  Barnes  (1973)  objective  analysis 
scheme.  Tests  were  conducted  to  determine  the  limitations 
of  the  Barnes  analysis  method  and  evaluate  how  well  this 
technique  reconstructs  meteorological  fields.  A  method  was 
developed  which  enables  an  optimum  weight  parameter  to  be 
determined  for  a  given  data  distribution. 

Barnes  objective  analyses  were  generally  more  accurate 
than  kinematic  fields  based  on  detailed  hand  analyses.  The 
objective  analysis  scheme  effectively  filtered  random  errors 
typical  of  instrumental  and  exposure  effects  from  the 
analysis.  However,  interpolation  near  the  edge  of  a  data 
domain  and  extrapolation  into  data  void  regions  led  to 


erroneous  results 


2 

An  optimum  weight  parameter  of  500  km  was  obtained 
in  a  mid-west  cyclone  test  case.  This  resulted  in  a  total 
response  of  0.936  for  wavelengths  of  twice  the  station 
spacing. 

The  Barnes  objective  analysis  method  was  also  used  in 
conjunction  with  mathematically  simulated  PAM  II  data  to 
estimate  the  probable  analysis  accuracy  of  the  system 
(analysis  plus  network).  Analyses  based  on  the  present 
PAM  II  network  design  did  not  adequately  represent  the  true 

meteorological  field  containing  mesoscale  structure. 

1 

Percentage  errors  exceeded  100%  in  the  case  of  temperature 
advection  and  local  frontogenesis.  However ^  the  proposed 
PAM  II  network  provided  analyses  nearly  as  good  as  those 
using  a  twice  as  dense  network,  and  somewhat  better  than 
those  obtained  using  the  regular  hourly  reporting  stations. 
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HARMS,  DEWEY  ELVIN.  Application  of  an  Objective  Analysis 
Scheme  to  Mesoscale  Observational  Network  Design.  (Under 
the  direction  of  DAVID  BARBER  and  GERALD  WATSON.) 

This  research  program  was  initiated  with  the  overall 
objective  of  estimating  the  probable  effectiveness  of  the 
proposed  surface  observational  network  (PAM  II)  which  will 
be  deployed  as  part  of  the  winter  1986  GALE  field  project. 
This  network  is  to  contribute  to  the  study  of  mesoscale 
features  within  developing  east  coast  cyclones,  and 
especially  the  well-known  coastal  front. 

To  accomplish  this  objective,  the  process  of  spatial 
weighted-averaging  interpolation  of  surface  data  was 
investigated  by  using  the  Barnes  (1973)  objective  analysis 
scheme.  Tests  were  conducted  to  determine  the  limitations 
of  the  Barnes  analysis  method  and  evaluate  how  well  this 
technique  reconstructs  meteorological  fields.  A  method  was 
developed  which  enables  an  optimum  weight  parameter  to  be 
determined  for  a  given  data  distribution. 

Barnes  objective  analyses  were  generally  more  accurate 
than  kinematic  fields  based  on  detailed  hand  analyses.  The 
objective  analysis  scheme  effectively  filtered  random  errors 
typical  of  instrumental  and  exposure  effects  from  the 
analysis.  However,  interpolation  near  the  edge  of  a  data 
domain  and  extrapolation  into  data  void  regions  led  to 


erroneous  results. 
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An  optimum  weight  parameter  of  500  km  was  obtained 
in  a  mid-west  cyclone  test  case.  This  resulted  in  a  total 
response  of  0.936  for  wavelengths  of  twice  the  station 
spacing . 

The  Barnes  objective  analysis  method  was  also  used  in 
conjunction  with  mathematically  simulated  PAM  II  data  to 
estimate  the  probable  analysis  accuracy  of  the  system 
(analysis  plus  network) .  Analyses  based  on  the  present 
PAM  II  network  design  did  not  adequately  represent  the  true 
meteorological  field  containing  mesoscale  structure. 
Percentage  errors  exceeded  100%  in  the  case  of  temperature 
advection  and  local  f rontogenesis.  However,  the  proposed 
PAM  II  network  provided  analyses  nearly  as  good  as  those 
using  a  twice  as  dense  network,  and  somewhat  better  than 
those  obtained  using  the  regular  hourly  reporting  stations. 
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CHAPTER  1 


INTRODUCTION 

"Objective  analysis"  is  a  programmable  method  for 
estimating  grid  point  values  by  transforming  (interpolating) 
data  from  observations  at  irregularly  spaced  points  into 
data  at  the  points  of  a  regular  grid  mesh.  The  gridded 
information  can  be  used  for  diagnostic  studies  or  in 
numerical  weather  prediction  models. 

The  method  of  analysis  is  one  of  the  the  several  factors 
which  determines  the  analysis  accuracy  achievable  with  a 
given  observational  network.  The  accuracy  also  depends  on 
the  observational  configuration,  density,  and  frequency;  the 
element  being  analyzed;  and  the  characteristics  of  the 
otiser  va  t  i  ona  1  error  (Seaman,  1977).  This  study  deals 
f  r  1  rr  a  r  1  1  y  with  the  method  of  analysis  and  with  the  density 
’  nfjquration  of  an  observational  network  to  investigate 
.  .  thp  tystem  (analysis  plus  network)  represents  true 

'  r  .  ;  ;  a  1  fields. 

,  r*'Sfarrh  was  conducted  in  support  of  the  East  Coast 
>  '■  r  le-t  (GALP)  to  estimate  the  probable 

•  r,f  the  proposed  surface  observational  network. 
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the  Portable  Automated  Mesonet  (PAM  II).  Three  of  the  basic 
objectives  of  GALE  are; 

--To  observe  and  document  the  structure  and  evolution  of 
mesoscale  weather-producing  circulation  systems  associated 
with  mid-Atlantic  coastal  cyclone  development; 

--To  quantify  the  physical  processes  which  lead  to 
development  of  these  circulation  systems,  in  particular, 
coastal  f rontogenesis ; 

--To  study  the  relationships  between  the  circulation 
systems  and  the  significant  weather  with  which  they  are 
associated . 

The  accomplishment  of  these  goals  will  require  the 
utilization  of  a  mesoscale  observing  network  which  describes 
the  circulation  systems,  that  is,  resolves  the  intense 
gradients  which  often  accompany  these  systems,  and  the 
application  of  improved  diagnostive  techniques  to  the 
comprehensive  data  sets  which  will  be  collected  by  this 
network  during  GALE. 

The  purpose  of  the  present  research  is  to  investigate 
the  process  of  spatial  weighted-averaging  interpolation  of 
surface  data  using  the  Barnes  (1973)  objective  analysis 
scheme.  Error  analyses  are  conducted  to  assess  resolution 
limits  of  an  observing  network,  in  particular,  the  PAM  II 
observation  system. 

The  combination  of  the  PAM  II  network  with  this 
objective  analysis  scheme  may  be  considered  as  a  tool  for 
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that  forces  a  high  degree  of  convergence  (agreement)  between 

the  observed  field,  f(x,y),  and  the  correction-pass 

interpolated  field,  g^(x,y).  Therefore,  and  y  can 

be  used  to  achieve  the  desired  response  of  waves  in  the 

analysis  (Koch  et  al.,  1981,  1983). 

With  several  tests  on  different  data  distributions, 

Barnes  (1973)  and  Koch  et  al.  (1981)  have  shown  that  >  of 

approximately  0.3  yields  the  best  results.  Specifying 

y  >  0.4  does  not  produce  rapid  analysis  convergence; 

conversely,  using  y  <  0.2  leads  to  excessive  "noise"  in  the 

analyzed  field.  Therefore,  a  value  of  0.3  is  used  in  all 

tests  described  in  this  thesis. 

For  computational  expediency  the  author  modified  the 

weight  function  W^.  The  weight  is  assigned  a  value  of 

zero  if  the  distance  between  the  grid  point  and  the  datum 

exceeds  a  chosen  critical  or  "influence"  radius  r  .  The 

c 

influence  radius  is  assigned  a  value  large  enough  (6-7  times 

as  large  as  the  average  station  spacing.  An)  which  insures 

that  only  data  which  would  otherwise  be  weighted  extremely 
—  8  0 

small  (W  <10  )  receive  zero  weight, 

n 

The  final  or  "correction-pass"  interpolated  grid  field, 
gj^(i,j),  is  the  sum  of  the  first  pass  field  and  the 
smoothed  residual  difference  between  the  observation  values, 
f(x,y),  and  the  first  pass  estimated  values  at  the  data 
location,  gQ(x,y).  A  simple  bilinear  interpolation 
between  the  values  of  g^(i,j)  at  the  four  nearest  grid 
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Figure  2.1  -  Relationship  of  first  pass  response  D  to  wave¬ 
length  ^  (km)  for  various  values  of  weight 
parameter  (from  Figure  1  of  Koch  et.  al  , 
1981).  ° 


•'o"  denotes  the  initial  pass.)  The  interpolated  field  at 
grid  points  is  thus  given  by 

=  f(x,y)D^  ,  (2.3) 

where  f(x,y)  is  an  idealized  observed  field,  that  is,  a 
continuum  field  with  no  errors.  It  is  important  to  note 
that  the  filter  response  function  that  Barnes  derived  is 
based  upon  an  idealized  field,  but  in  reality  only  finite 
sets  of  observed  data  with  some  degree  of  random  error  are 
available.  Fig.  2.1  reveals  that  the  incorporation  of 
smaller  values  of  the  first  pass  weight  parameter,  k^, 
results  in  greater  filter  response,  particularly  for  the 
short  waves. 

The  original  (1964)  version  of  the  interpolation  scheme 
required  several  iterations  through  the  data  field.  Also, 
the  weight  parameter  was  the  same  for  all  passes.  Barnes 
(1973)  modified  the  earlier  (1964)  version  of  his  technique 
in  order  to  decrease  the  computation  time  necessary  to 
achieve  the  desired  response  at  small  wavelengths.  This 
modification  enables  a  single  correction  pass  through  the 
initial  interpolated  field  to  achieve  the  desired  response 
for  specific  wavelengths,  rather  than  requiring  several 
iterations  as  before.  This  is  accomplished  by  decreasing 
the  weight  parameter,  k,  from  its  first  pass  value  of  k^ 
to  a  second  (correction)  pass  value  of 

k,  =  >k^  .  (2.4) 

Here  >  is  a  "numerical  convergence  parameter",  0  <  "r  <  1 , 
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indicate  that  the  Barnes  analysis  technique  is  as  good  as 
detailed  hand  analyses.  The  inter-comparisons  reveal  that 
the  differences  among  the  analyses  are  of  the  same  order  as 
those  when  the  comparison  is  made  with  the  true  data  set. 

2.3  Characteristics  of  the  Barnes  Objective  Analysis  Method 

Barnes  (1964)  developed  an  objective  analysis  method 

that  uses  a  Gaussian  weight  function  in  the  spatial  domain 

based  on  the  fundamental  premise  that  two-dimensional  data 

distributions  can  be  represented  by  Fourier  integrals. 

Weights  are  assigned  to  the  data  as  a  known  function  of 

distance  between  observations  and  a  particular  grid  point. 

The  weight,  is  assigned  according  the  distance, 

d^  (in  kilometers,  km),  between  the  grid  point  (i,j)  and 

the  datum  f  (x,y)  as 
n 

=  exp(-d^/4k)  ,  (2.1) 

where  k  is  a  weight  parameter  that  determines  the  shape  of 

the  filter  response  function  (the  percentage  of  the  original 

wave's  amplitude  that  will  be  conserved).  (NOTE:  k  is  in 
2 

units  of  km  which  will  be  the  assumed  units  for  weight 
parameters  throughout  this  thesis.) 

With  one  pass  through  the  data,  the  filter  response 
function  can  be  shown  to  be 

=  exp  (-k^  (21' /X  )  ^)  ,  (2.2) 

where  X  is  the  the  horizontal  wavelength.  (The  subscript 

» 
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Table  2.1  -  Results  of  a  comparison  of  the  Barnes  analysis 
with  hand  analyses  using  simulated  meteorologi¬ 
cal  fields.  Root-roean-square  differences  among 
various  analyses  are  presented. 


Root-Mean-Square  Difference 
Surface  Temperature  1000  mb  Height 


Barnes  /  Actual 

0.97  °F 

2.38  m 

Analyst  #1 

/ 

Actual 

0.79 

2.67 

Analyst  #2 

/ 

Actual 

0.43 

3.22 

.jnalyst  #1 

/ 

Barnes 

0.91 

3.33 

Analyst  #2 

/ 

Barnes 

1.11 

3.61 

Analyst  #1 

/ 

Analyst  #2 

0.78 

2.20 

pattern  scales  resolvable  by  the  data  distribution  will  be 
revealed . 

3. )  Since  the  weighting  function  approaches  zero 
asymptotically,  the  influence  of  data  may  be  extended  any 
distance  to  insure  that  a  sufficient  number  of  observations 
influence  each  grid  point  without  changing  the  weight 
function  and,  therefore,  the  response  characteristics. 

4. )  Only  two  passes  through  the  data  are  needed  to 
achieve  the  desired  pattern. 

5. )  "Noise"  is  adequately  filtered  from  the  analysis  so 
further  smoothing  by  application  of  additional  numerical 
filters  is  not  necessary. 

In  support  of  the  choice  of  the  Barnes  analysis 
technique,  tests  were  conducted  to  compare  this  analysis 
method  with  subjective  hand  analyses.  The  hand  analyses 
were  performed  by  two  meteorology  graduate  students  with 
analysis  experience.  "Bogus"  or  simulated  fields  of  surface 
temperature  and  1000  mb  height  were  analyzed.  (Construction 
of  these  fields  through  utilization  of  mathematical 
functions  will  be  discussed  at  length  in  Chapter  4.)  The 
resulting  grid  point  values  from  the  two  hand  analyses  and 
the  Barnes  scheme  were  compared  to  the  "true"  values  at 
these  locations,  which  were  calculated  directly  with  the 
simulation  functions.  Inter-comparisons  were  also  made 
among  the  three  analyses.  The  results  of  the  comparisons, 
shown  as  root-mean-square  differences  (RMSD)  in  Table  2.1, 


shape  of  the  surface. 

Two  of  the  most  used  successive  correction 
(weighted-averaging)  techniques  are  the  Barnes  (1973)  and 
Cressman  (1959)  objective  analysis  schemes.  At  first 
glance,  the  Cressman  scheme  may  seem  more  attractive  since 
it  does  not  incorporate  a  Gaussian  weight  function  and  is 
therefore  especially  cost  effective.  However,  in  this 
technique,  the  weight  function  is  dependent  on  the  influence 
radius  beyond  which  zero  weight  applies;  the  weights  do  not 
approach  zero  asymptotically  as  in  the  Barnes  scheme.  To 
insure  that  sufficient  data  influence  the  interpolation  in 
data  sparse  regions,  the  current  scan  radius  is  increased 
locally.  This  procedure  produces  un)cnown  response  in  the 
final  result  and  introduces  short  wave  features  (noise) 
which  must  be  smoothed  by  later  application  of  arbitrary 
filters.  Additionally,  four  or  more  passes  through  the  data 
are  required  in  order  that  Cressman 's  method  achieve  the 
desired  pattern. 

To  summarize,  the  Barnes  (1973)  scheme  is  selected  for 
objective  analyses  in  this  study  because  of  the  following 
advantages  (listed  by  Barnes) : 

1. )  The  scheme  is  computationally  simple  and  assigns 
weights  solely  as  a  function  of  distance  between  datum  and 
grid  point. 

2. )  The  weight  parameter,  which  determines  the  degree  of 
smoothing,  can  be  chosen  prior  to  the  analysis  so  that 
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in  this  study  because  of  its  simplicity,  accuracy,  cost 
effectiveness,  and  its  several  advantages  (mentioned  later) 
over  similiar  techniques.  The  method  of  statistical 
interpolation  was  not  chosen  since  this  technique  is 
complex,  and  also  tends  to  exhaust  computer  resources  when 
calculating  the  observational  and  initial  guess  correlations 
which  determine  the  weighted  coefficients.  Furthermore, 
Otto-Bliesner  et  al.  (1977)  showed  that  this  method  is  no 
better  or  worse  than  methods  similar  to  the  Barnes  scheme 
(successive  corrections)  when  the  analyses  are  performed 
over  a  data-rich  region.  The  variational  technique  of 
Sasa)ci  (1958,  1970)  also  attempts  to  accomplish  not  only  the 
interpolation  phase  of  objective  analysis  but  also  the 
matching  of  variables  (data  assimilation)  which  requires 
additional  computation  time.  For  the  purposes  of  this 
research,  only  the  interpolation  of  observational  data  to  a 
regular  grid  is  necessary. 

In  the  past  the  method  of  polynomial  (surface)  fitting 
has  been  successful  in  accomplishing  the  desired 
interpolation.  However,  according  to  Barnes  (1964),  such 
schemes  possess  three  major  disadvantages:  the  calculations 
are  complicated  and  costly  with  respect  to  computer  time; 
the  data  to  which  the  surface  is  fitted  are  chosen  in  a 
rather  artificial  manner  (that  which  produces  the  best 
results) ;  and  the  effect  of  erroneous  data  can  be  disastrous 
since  each  datum  is  given  equal  ranking  in  determining  the 
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minimum;  also,  they  depend!  upon  the  spatial  covariances 
among  the  analyzed  variables.  Examples  of  this  technique 
are  found  in  Gandin  (1963),  Eddy  (1967),  Bengtsson  and 
Gustavsson  (1971),  Rutherford  (1972,  1973),  Thiebaux  (1975), 
and  Schlatter  et  al.  (1976). 

3. )  Variational  Technique.  The  method,  developed  by 
Sasalci  (1958,  1970),  is  a  post-analysis  adjustment  technique 
based  upon  the  calculus  of  variations.  A  functional  is 
utilized  which  minimizes  the  analyzed-minus-observed 
differences,  filters  undesirable  high-frequency  and 
high-wave  number  features,  and  employs  dynamical 
constraints.  These  constraints  may  be  strong  (satisfied 
exactly)  or  wea)c  (satisfied  only  approximately).  This 
method  is  effective  in  ma)ting  the  analyzed  fields  compatible 
with  a  forecast  model, 

4. )  Polynomial  Fitting.  Observations  are  fitted  by 
utilizing  polynomials  of  two  or  more  variables.  The  fit  is 
obtained  by  a  least  squares  technique  since  there  are  often 
fewer  polynomial  coefficients  than  data  points.  Polynomial 
fitting  techniques  have  been  used  by  Panofs)cy  (1949), 
Gilchrist  and  Cressman  (1954),  Penn  et  al.  (1963),  and 
Shapiro  and  Hastings  (1973). 

2.2  Comparison  of  the  Barnes  and  Alternative  Schemes 


The  Barnes  objective  analysis  scheme  was  chosen  for  use 


To  review  the  different  methods  of  objective  analysis, 
four  of  the  more  prominent  techniques  are  summarized 
(Schlatter  et  al . ,  1976): 

1. )  Successive  Corrections.  A  forecast  model  or  an 
interpolation  from  first  pass-grid  values  back  to  station 
locations  provides  the  preliminary  estimate  (first-guess)  of 
the  field  to  be  analyzed.  This  estimate  is  modified  by  a 
combination  of  corrections  computed  for  each  grid  point. 

The  corrections  are  proportional  to  the  differences  between 
the  observed  and  first-guess  values.  These  are  weighted 
empirically,  with  observations  nearest  the  grid  point 
weighted  most  heavily.  Several  scans  or  passes  are  made 
through  the  data;  the  results  of  one  scan  become  the  first 
guess  for  the  next  scan.  Bergthdrsson  and  Doos  (1955), 
Cressman  (1959),  and  Barnes  (1964,  1973)  have  developed 
analysis  schemes  based  upon  the  successive  corrections 
method . 

2. )  Statistical  (Optimum)  Interpolation.  This  analysis 
technique,  based  on  statistical  linear  regression, 
incorporates  either  a  forecast  model,  persistence,  or 
climatology  as  a  first-guess  field.  The  analyzed  value  at  a 
grid  point  is  the  sum  of  this  first-guess  and  a  linear 
combination  of  corrections,  which  are  proportional  to  the 
differences  between  observational  and  first-guess  values. 

The  weighting  coefficients  are  determined  from  the  condition 
that  the  mean-square-error  of  the  analyzed  values  be  a 


smoothing  is  essential  to  producing  fields  that  will  provide 
meaningful  estimates  of  such  derived  quantities  as 
divergence,  advections,  and  gradients  from  the  grid  point 
values  (Koch  et  al.,  1981). 

As  noted  by  Gandin  (1963),  the  term  "objective  analysis" 
is,  to  some  extent,  misleading.  The  analysis  of  a 
meteorological  field  does  not  imply  the  reconstruction  of 
the  whole  field.  Actually  it  implies  a  study  of  the 
properties  of  an  already  constructed  field.  Thus,  according 
to  Gandin,  "numerical  methods  for  the  reconstruction  of 
meteorological  fields"  would  be  a  more  appropriate  term. 
However,  the  term  "objective  analysis"  will  be  used  in  this 
text  since  it  is  widely  accepted  and  used. 

Objective  analysis  performs  three  basic  tas)cs:  1.)  the 
interpolation  of  values  of  observed  meteorological  elements 
to  grid  points,  2.)  the  assimilation  of  raw  data,  a  process 
in  which  meteorological  observations  are  absorbed  into  a 
numerical  representation  of  the  atmosphere,  and  3.)  the 
detection  and  elimination  of  erroneous  data.  "Assimilation 
may  be  distinguished  from  interpolation  by  requiring  that 
the  numerical  representation  be  a  complete  description  of 
meteorological  variables  whose  mutual  dependence  is  governed 
by  explicit  physical  constraints"  (McPherson,  1975).  Of  the 
three  basic  tasks  listed,  the  first  (interpolation)  is  the 
most  important  part  of  objective  analysis  and  will  be 
investigated  at  length  in  this  paper. 


5 


CHAPTER  2 

BARNES  OBJECTIVE  ANALYSIS  METHOD 

Before  an  in-depth  description  of  the  Barnes  objective 
analysis  method  is  presented,  it  is  appropriate  to  first 
review  the  purpose  and  content  of  objective  analysis,  and  to 
provide  a  brief  summary  of  the  better  known  methods. 

2.1  Review  of  Objective  Analysis  Methods 

Objective  analysis  involves  the  development  of 
mathematical  methods  which  use  data  at  irregularly  spaced 
points  to  reconstruct  objectively  the  continuous  field  of 
some  meteorological  variable,  and  to  display  this  field  at  a 
regular  set  of  grid  points.  The  data  in  this  format  then 
permit  quantitative  diagnostic  and/or  predictive 
calculations  to  be  made. 

The  process  of  objective  analysis  results  in  the 
smoothing  (suppressing)  of,  especially,  the  higher 
frequency,  shorter  wavelength  features,  including 
unavoidable  "noise”,  and  those  arising  from  large  data 
errors  or  nonrepresentative  observations.  Controlled 
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through  the  utilization  of  simulated  meteorological  fields, 
the  influence  of  this  analysis  technique  on  the 
reconstructed  fields.  In  addition,  an  estimate  is  made  of 
the  confidence  with  which  kinematic  fields  can  be  evaluated 
by  superposing  random  errors  typical  of  instrumental  and 
site  exposure  variability  on  the  simulated  field.  This 
procedure,  to  some  extent,  evaluates  the  limitation  of  a 
meteorological  network  in  measuring  or  describing  mesoscale 
weather  systems.  Furthermore,  these  tests  can  be  helpful  in 
making  final  adjustments  to  the  proposed  PAM  II  network 
design,  and  possibly  in  developing  strategies  for  the 
deployment  of  mobile  observing  systems  for  the  GALE  project. 
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obtaining  estimates  of  meteorologically  significant 
parameters  such  as  vorticity,  divergence,  gradients,  and 
advections.  To  obtain  an  estimate  of  how  well  that  tool 
will  function,  the  Barnes  objective  analysis  scheme  is 
applied  to  bogus  data  sets  which  include  a  reasonable  error 
component.  These  data  sets  are  developed  by  judicious 
enhancement  of  case  study  results  using  a  simulation  model 
of  known  mathematical  form,  and  dependent  only  on  a 
station's  location.  Using  existing  studies  of  coastal 
fronts,  gradients  are  sharpened  beyond  those  which  can  be 
resolved  by  the  observation  network  while  not  violating  any 
of  the  data.  Values  (from  the  resulting  fields)  which 
correspond  to  the  proposed  PAM  II  station  locations  are  used 
as  simulated  observational  data.  These  "observations"  based 
on  the  known  field  are  then  made  more  realistic  by  adding  a 
reasonable  error  component.  Finally,  these  data  are 
analyzed  to  a  mesoscale  grid.  The  resulting  derived  fields 
are  compared  with  those  obtained  without  the  introduction  of 
error  and  to  the  "true"  simulated  field  values  at  the 
corresponding  grid  points. 

More  specifically,  the  task  of  this  research  effort  is 
to  utilize  a  highly  flexible  version  of  the  Barnes  objective 
analysis  scheme  to  perform  tests  to  determine  optimum  weight 
parameters  (which  control  the  degree  of  smoothing  performed 
in  the  analysis)  to  be  used  in  the  study  of  mesoscale 
processes.  An  attempt  is  also  made  to  isolate  and  study. 


points  is  used  to  estimate  g^(x,y). 
given  by 


The  final  field  is 
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g],(i»j)  =  (i»  j)  +  lf  (x,y)-g^  (x,y)  ]Dj^  .  (2.5) 


Here 


Dj  =  exp[-kj^  (2ti/X  )  ^] 

=  exp[-Yk^(2r/A)^J  =  (2.6) 

is  the  response  function  corresponding  to  the  weight 
function 

”n  '  exp(-d2/4.k„)  ,  a„  < 

0  *  ^  '^c  ‘ 

The  actual  correction-pass  value  at  each  grid  point  results 
from  adding  the  weighted  averages  from  the  two  passes  with  N 
observations  (This  follows  from  2.5.): 


g^  (i» j) 


E 

T^=  1 


N 


n=l 


n=l  n=l 


(2.8) 


Substituting  (2.6)  and  (2.3)  into  (2.5)  gives 
gj^(i,j)  =  D^f  (x,y)  +  (f  (x,y)-D^f  (x,y)  ]Dq 

=  f  (x,y)D^(l  +  D^‘^-D^)  ;  (2.9) 

therefore,  the  true  correction  pass  (total)  response 
function  is 

This  total  response  function,  D^,  indicates  the  degree 
of  analysis  convergence,  that  is,  how  closely  the 
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interpolated  values  agree  with  the  observed  ones.  Fig.  2.2 
shows  that  the  second  pass  through  the  data  will  increase 
the  degree  of  convergence  when  0  <  >  <  1. 

Koch  et  al .  (1981)  have  proven  that  the  1973  version  of 
the  Barnes  objective  analysis  technique  is  absolutely 
convergent.  This  fact  enables  the  analyst  to  control  the 
amount  of  small-scale  detail  to  be  revealed  in  the  analyzed 
data  fields. 

2.4  Computer  Program  Package 

The  author  has  developed  a  computer  program  which 
instills  versatility  in  the  objective  analysis  through  human 
selection  of  input  parameter  values:  grid  spacing  and 
origin,  weight  parameter,  and  numerical  convergence 
parameter.  Before  the  analysis  computations  are  performed, 
a  rectangular  coordinate  system  with  an  array  of  grid  points 
is  constructed  using  a  polar  stereographic  map  projection 
(Inman,  1970).  The  grid  can  be  constructed  over  any  area  of 
the  earth  at  any  desired  orientation.  A  highly  flexible 
input  routine  allows  easy  manipulation  of  raw  data.  The 
objective  analysis  procedure  follows  the  logic  outlined  in 


Section  2.3. 


20 


2.5  Research  Goals 

Through  the  use  of  the  Barnes  objective  analysis  scheme, 
several  goals  of  this  research  are  achieved.  The  goals 
include : 

1. )  Determination  of  "edge”  effects;  that  is,  how  the 
analysis  scheme  behaves  near  the  edge  of,  or  beyond  the  data 
domain . 

2. )  Determination  of  the  optimum  weight  parameter  and 
grid  spacing  for  a  given  data  distribution. 

3. )  Determination  of  the  effects  of  incorporating  random 
error  in  a  data  field. 

4. )  Determination  of  the  Barnes  scheme's  effects  on  the 
character  of  the  information  it  processes. 

5. )  Obtaining  an  a  priori  estimate  of  typical  errors 
which  may  occur  in  the  basis  and  derived  kinematic  fields 
based  on  GALE  PAM  II  data. 

6. )  Determination  of  whether  the  proposed  PAM  II  network 
produces  analyses  significantly  better  than  those  obtained 
from  the  regular  hourly  reporting  station  network. 

7. )  Experimentation  with  alternative  PAM  II  networks  to 
test  their  impact  on  the  depiction  of  features  of  particular 
interest  to  GALE  project  goals. 

The  following  chapter  describes  a  test  case  in  which  the 
first  three  objectives  are  investigated. 
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CHAPTER  3 

EXPERIMENTS  WITH  THE  BARNES  ANALYSIS 
METHOD  ON  A  MID-WEST  CYCLONE 

3.1  Overview  and  Synoptic  Situation 

The  case  of  9  November  1983  (2100  GMT)  over  the  central 
United  States  was  used  to  test  the  Barnes  objective  analysis 
computer  program,  and  to  investigate  "edge”  effects,  optimum 
weight  parameter,  effects  of  random  error  on  data  fields, 
and  overall  accuracy  as  compared  to  subjective  analyses. 

The  mid-west  test  case  provided  non-uniform  fields  and 
strong  gradients  needed  to  test  the  limitations  of  the 
Barnes  objective  analysis  and  determine  how  well  this 
technique  reproduces  true  meteorological  fields. 

This  particular  case  features  a  mature  surface  cyclone 
with  a  sharply-defined  frontal  system  over  the  Central 
Plains  (Fig.  3.1).  To  be  noted  are  the  pronounced 
temperature  gradient  and  cyclonic  shear  through  the  frontal 
zone  and  the  strong  pressure  gradient  west  of  the  cyclone 
with  a  rather  weak  gradient  to  the  east.  Since 
precipitation  was  occurring  at  the  chosen  synoptic  time 
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Figure  3.1a  -  Hand  analysis  of  surface  temperature  (°F)  for 
2100  GMT  9  November  1983.  Contours  are  drawn 
with  interval  of  4  °F. 


Figure  3.1b  -  Hand  analysis  of  sea  level  pressure  (mb)  for 
2100  GMT  9  November  1983. 
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Figure  3.1c  -  Hand 
wind 
Wind 
with 


plot  of  a  representative  sample  of  the 
reports  for  2100  GMT  9  November  1983. 
barbs  are  plotted  using  standard  notation 
speed  in  knots. 


(2100  GMT),  a  qualitative  comparison  of  these  areas  can  be 
made  with  the  objectively  analyzed  patterns  of  surface 
convergence  and  relative  vorticity. 

The  objective  analysis  was  performed  over  a  rectangular 
area  whose  sides  extend  from  North  Dakota  to  the  Texas 
Panhandle  to  northern  Georgia  to  the  upper  Great  Lakes 
region.  A  total  of  143  station  observations  were 
incorporated  into  the  analysis.  A  15X18  grid  was  constructed 
on  a  regular  coordinate  system  with  the  x-  axis  along  the 
central  longitude  of  the  grid  positive  to  the  south  and  the 
y-  axis  positive  to  the  east  at  the  center  of  the  grid.  A 
grid  spacing  of  115  km  was  used.  An  inner  grid  encompassing 
the  central  35  stations  of  the  total  network  was  also  used 
in  various  tests. 

A  qualitative  comparison  of  Barnes  and  hand  analyses  of 
pressure  and  temperature  patterns  can  be  made  by  examining 
Figs.  3.1  and  3.2.  Fig.  3.2c  shows  the  wind  field  based  on 
an  analysis  of  the  wind  components.  (The  u-  and  v- 
components  of  the  observed  winds  were  analyzed  seperately; 
then  these  objectively  analyzed  grid  values  of  u-  and  v- 
were  added  to  obtain  wind  speeds  and  directions  at  the  grid 
points.)  The  objectively  analyzed  wind  field  may  be 
compared  with  the  plotted  winds  in  Fig.  3.1c. 

Fig.  3.3a  shows  the  precipitation  reports  superimposed 
on  the  surface  divergence  pattern  produced  from  the  Barnes 
analyzed  u-  and  v-  wind  components.  The  precipitation 
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Figure  3.2a  -  Barnes  analysis  of  surface  temperature  (°F) 
for  2100  GMT  9  November  1983  (k^  =  500). 


Figure  3.2b  -  Barnes  analysis  of  sea  level  pressure  (mb) 
for  2100  GMT  9  November  1983  (k  =  500). 


Fiqure  3.2c  -  Vector  display  of  Barnes  analysis  of  surface 
winds  for  2100  GMT  9  November  1983  (k  =  500 
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Figure  3.3a  -  Surface  divergence  (10  /s)  derived  from 

Barnes  analyzed  wind  components  for  2100  GMT 
9  November  1983.  Precipitation  reports  are 
super  imposed . 


Figure  3.3b  -  Vertical  component  of  relative  vorticity 

(10~^/s)  derived  from  Barnes  analyzed  wind 
components  for  2100  GMT  9  November  1983. 
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pattern  approximately  coincides  with  the  axis  of  maximum 
convergence.  The  objective  analysis  of  surface  relative 
vorticity  (Fig.  3.3b)  shows  maximum  positive  vorticity  along 
the  trough  and  at  the  low  pressure  center,  as  one  would 
expect . 

In  each  of  the  above  figures  the  graphic  package 
(SURFACE  II),  developed  at  the  University  of  Kansas,  was 
used  to  objectively  contour  the  gridded  data  (Figs.  3.2  and 
3.3).  Thus,  the  analyses  are  produced  with  complete 
objectivity  by  the  computer. 

3.2  Edge-Effect  Test 

It  would  be  advantageous  to  an  analyst  to  know  how  the 
analysis  scheme  behaves  near  the  edge  of  the  data  field. 

With  this  information  one  can  decide  whether  to  retain  or 
reject  the  gridded  values  near  the  edge  of  or  beyond  the 
data  domain.  If  these  values  are  accepted,  one  will  have  a 
better  insight  into  their  reliability.  The  following  steps 
were  taken  to  help  determine  these  effects: 

1. )  The  Barnes  analysis  was  performed  over  the  entire 
grid  using  all  143  stations.  This  analysis  was  used  to 
provide  standard  or  "true"  values  over  the  entire  region. 

2. )  The  objective  analysis  was  then  performed  over  an 
inner  array  (centered  within  the  larger  grid)  of  11X11  grid 
points  and  using  35  stations.  An  inner-most  7X7  grid  array 
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approximately  coincided  with  the  region  encompassing  these 
stations. 

3.)  Each  of  the  grid  point  values  from  the  smaller  grid 
was  subtracted  from  the  corresponding  grid  value  in  the 
larger  grid.  Root-mean-square  errors  (RMSE)  were  calculated 
using  the  values  of  the  peripheral  grid  points  of  each 
successive  smaller  square  of  grid  points,  where  N  =  40,  32, 
24,  16,  and  8  for  the  successive  squares,  using: 

RMSE  =  (GPj^-GPg)  (3.1) 

where  GP,  is  the  grid  point  value  from  the  15X18  grid 

L« 

and  GP  is  the  corresponding  value  from  the  inner  grid. 

The  results,  shown  in  Table  3.1,  reveal  a  significant 
increase  in  error  near  the  edge  of  the  data  field  with 
approximately  a  200-300%  increase  in  error  from  the  edge  of 
the  data  domain  (7X7)  to  the  outer  square  (11X11).  Hence, 
it  is  shown  that  extrapolation  into  data  void  regions  leads 
to  extremely  erroneous  results  and  even  as  the  edge  of  the 
data  domain  is  approached  the  interpolation  becomes 
unrel iable . 

3.3  Weight  Parameter  and  Random  Error  Tests 

The  choice  of  a  weight  parameter  for  a  given  data 
distribution  is  critical  since  the  weight  determines  the 
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Table  3.1  -  "Edge  Effects  Test":  Root-mean- square-error  for 
successive  peripheral  squares  of  grid  points  for 
which  the  15X18  grid  serves  as  the  standard. 
Total  number  of  grid  points  is  given  for  each 
square.  Results  for  temperature  and  wind  com¬ 
ponents  are  shown. 


11X11  sq. 

9X9 

7X7 

5X5 

3X3 

(40  pts.) 

(32) 

(24) 

(16) 

(8) 

TEMP  (°F) 

2,80 

1.59 

1.06 

0.89 

0.63 

U  (m/s) 

3.71 

2.22 

1.10 

0.48 

0.26 

V  (m/s) 

1.03 

0.93 

0.76 

0.52 

0.30 
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CHAPTER  4 

A  SIMULATION  OF  METEOROLOGICAL  FIELDS 

4.1  Construction  of  the  Analytic  Functions 

To  determine  how  well  a  particular  arrangement  of  the 
PAM  II  network  will  provide  estimates  of  meteorologically 
significant  parameters  such  as  divergence,  vorticity, 
temperature  gradient  and  advection,  experimental  "bogus" 
data  sets  will  be  used  with  the  Barnes  objective  analysis. 
Meaningful  tests  of  this  analysis  method  require  data  sets 
at  both  grid  and  observation  points.  A  simulation  model 
incorporating  mathematical  functions  which  are  dependent  on 
a  point's  latitude  and  longitude  is  chosen  to  generate 
meteorological  fields  of  temperature  and  geopotential 
height.  These  analytical  fields  should  resemble  actual 
meteorological  fields  and  possess  enough  non-linearity  to 
identify  problems  and  merits  of  the  observing  network  and 
interpolation  scheme  being  tested. 

The  analytic  functions  chosen  are  based  upon 
Sanders'  (1971)  structural  model  of  temperature  and  height 
fields  which  he  used  to  find  analytic  solutions  to  the 


Table  3.6  -  Effects  on  optimuin  weight  parameter  and  RMSEs 
due  to  increasing  the  grid  spacing.  (Same  as 
Table  3.5b  except  DX  =  115  km.) 


(Divergence  and 
wind  components 

vorticity 
in  units 

in 

of 

units  of  10 
m/s) 

Vs, 

Root-Mean- 

Square-Error 

‘'o 

DIV 

VORT 

U 

V 

12.5 

0.93 

0.85 

1.44 

1.20 

25. 

0.83 

0.62 

1.01 

0.99 

50. 

0.80 

0.61 

0.85 

0.76 

100. 

0.80 

0.85 

0.98 

0.50 

250. 

1.09 

1.19 

1.22 

0.54 

500. 

1.51 

1.41 

1.70 

0.76 

1000. 

1.91 

1.80 

2.34 

0.92 

1500. 

2.16 

2.09 

2.75 

0.99 

2000. 

2.22 

2.30 

3.03 

1.03 

3000. 

2.35 

2.54 

3.39 

1.08 

5000. 

2.47 

2.75 

3.73 

1.13 

-3  2 

(DX  =  115  km)  twice  as  large  as  before  to  determine  if  there 
is  a  significant  difference  in  the  mean  errors  and  the 
optimum  weight  parameter  of  the  kinematic  fields  when  DX  is 
varied.  when  the  results  of  Table  3.6  (DX  =  115  km)  are 
compared  to  those  of  Table  3.5b  (DX  =  57  km),  it  is  readily 
apparent  that  the  optimum  weight  of  each  of  the  four 
variables  remained  the  same  except  in  the  case  of  vorticity, 
where  the  optimum  weight  parameter  decreased  slightly.  The 
minimum  error  did  not  increase  a  significant  amount  with  the 
coarser  grid;  in  fact,  the  error  in  vorticity  actually 
decreased  somewhat.  These  results  concur  with  the  findings 
of  the  test  in  Section  3.3.1  which  incorporated  the  "missing 
station"  procedure.  Therefore,  it  is  concluded  that  using  a 
grid  spacing  significantly  smaller  than  the  average  station 
spacing  (DX  <  An/2)  is  not  justified  when  computing  costs 


are  considered 


The  results  from  the  more  reliable  test  method  (more 


reliable  since  the  standard  was  the  actual  known  station 

value)  in  Section  3.1.1  showed  the  optimum  weight  parameters 

for  the  u-  and  v-  wind  components  to  be  approximately  250 

and  500,  respectively  (Table  3.2).  This  suggests  that  the 

optimum  weight  parameters  determined  using  hand  analyses  as 

the  standard  may  be  biased  toward  too  small  values. 

Therefore,  if  the  optimum  weight  for  divergence  and 

vorticity  is  again  taken  to  be  twice  as  large  as  that  for  u- 

and  approximately  equal  to  that  of  v-,  the  appropriate 

optimum  k  for  the  kinematic  variables  would  be 
o 

approximately  500  for  the  given  station  distribution. 

With  an  average  station  spacing  of  approximately  95  km 
and  of  500,  the  total  response  (Equation  2.10)  of  the 
2  An  wave  is  found  to  be  0.936.  With  this  total  response  an 
"optimum”  weight  parameter  can  be  calculated  (using 
Equations  2.2  and  2.10)  for  similar  data  distributions  with 
known  average  station  spacings. 

The  calculated  response  of  0.936  may  seem  rather  high, 
especially  compared  to  Koch's  (1981)  suggestion  that  this 
response  should  equal  1/e  or  0.3679.  However,  to  achieve  a 
total  response  of  1/e  for  2  An  waves  and  y  =  0.3,  k^ 
would  have  to  be  approximately  3200.  From  the  results  in 
Fig.  3.4  it  is  apparent  that  the  error  would  be  unacceptably 
large,  primarily  due  to  over- smoothi ng . 

This  test  was  repeated  using  a  grid  spacing 


•;o 

shows  the  RMS  differences  of  each  individual  analysis 
compared  to  the  average  of  the  other  four  analyses.  The 
average  of  all  five  subjective  analyses  served  as  the 
standard  for  comparison  with  the  objectively  analyzed 
var iables . 

The  observed  u-  and  v-  wind  component  fields  were 

objectively  analyzed  using  various  weight  parameters.  These 

analyzed  fields  were  then  used  to  compute  divergence  and 

vorticity  fields.  Table  3.5b  shows  the  error  for  each 

variable  at  all  k  *s  tested.  The  RMSE  of  the  optimum 

o 

weight  for  each  variable  was  almost  as  low  as  the  error  of 
the  best  hand  analysis  of  the  corresponding  variable, 
particularly  for  divergence  and  vorticity.  For  these 
derived  fields  the  minimum  RMSE  was  less  than  the  majority 
of  the  hand  analysis  results.  Unfortunately,  the  error  in 
the  divergence  and  vorticity  approaches  the  same  order  of 
magnitude  as  the  values  themselves. 

The  results  in  Table  3.5b  enable  some  conclusion  to  be 
drawn  on  the  choice  of  an  optimum  weight  parameter.  The 
optimum  for  divergence  and  vorticity  is  about  100  and 
is  approximately  equal  to  that  for  v-  and  approximately 
twice  that  for  the  u-  wind  component  (optimum  k^  =  50). 

This  suggests  that  the  original  fields  (temperature,  wind 
components,  etc.)  may  require  smoother  analyses  when 
kinematic  fields  ( f ini te-d i f f erenced  values)  are  to  be 


constructed . 
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Table  3.5a  -  Root-mean- squa re  differences  of  each  hand 

analysis  compared  to  the  average  of  the  other 
four  analyses  with  DX  =  57  km.  Results  for 
surface  divergence,  relative  vorticity,  and  wind 
components  are  shown. 

_  5 

(Divergence  and  vorticity  in  units  of  10  /s  , 
wind  components  in  units  of  m/s) 


DIV 

VORT 

U 

V 

Analyst 

#1 

1.44 

1.02 

0.86 

0.96 

Analyst 

#2 

0.98 

0.79 

0.60 

0.47 

Analyst 

#3 

1.10 

0.78 

0.73 

0.30 

Analyst 

#4 

0.76 

0.80 

0.76 

0.38 

Analyst 

#5 

1.48 

1.11 

0.75 

1.16 

Table  3.5b  -  Weight  parameter  test  using  kinematic  variables 
with  DX  =  57  km.  Root-mean-square-errors  for 
each  weight  parameter  are  presented  for  which 
the  average  of  the  five  hand  analyses  served  as 
the  standard. 


DIV 

VORT 

U 

V 

12.5 

1.99 

2.07 

1.31 

0.91 

25. 

1.49 

1.59 

0.90 

0.90 

50. 

1.03 

1.10 

0.79 

0.73 

100. 

0.77 

0.78 

0.97 

0.48 

250. 

1.03 

1.14 

1.23 

0.54 

500. 

1.47 

1.39 

1.71 

0.76 

1000. 

1.91 

1.81 

2.34 

0.92 

1500. 

2.12 

2.13 

2.75 

0.99 

2000. 

2.25 

2.35 

3.03 

1 . 03 

3000. 

2.39 

2.62 

3.38 

1.07 

5000. 

2.52 

2.86 

3.73 

1.12 
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The  comparison  of  Table  3.4  (with  error)  and  Table  3.3 

(without  error)  reveals  the  optimum  weight  parameter  for 

each  variable  remains  unchanged.  Unexpectedly,  the  RMSE 

corresponding  to  the  optimum  decreases  slightly  when 

random  errors  are  added  to  the  temperature  and  v-  fields. 

However,  with  several  of  the  non-optimum  k  's  the  error 

o 

increases,  as  would  be  expected.  in  any  event,  the  results 
for  all  the  variables  indicate  that  the  Barnes  objective 
analysis  effectively  eliminates  the  effects  of  random 
errors. 

3.3,3  Tests  employing  kinematic  variables 

A  weight  parameter  should  be  chosen  which  minimizes  the 
error  in  the  kinematic  fields  such  as  temperature  advection, 
divergence,  vorticity,  and  horizontal  moisture  convergence 
since  these  are  especially  meteorologically  significant; 
that  is,  they  are  directly  related  to  the  genesis  and 
maintenance  of  weather  systems.  Therefore,  a  test  to 
determine  the  optimum  weight  parameter  for  these  computed 
quantities  is  presented  here. 

Analyses  of  the  u-  and  v-  wind  components  performed  by 
five  NCSU  faculty  and  graduate  students  served  to  produce  a 
standard  against  which  to  compare  the  Barnes  analysis 
results.  Divergence  and  vorticity  values  were  calculated 
using  these  analyses  of  the  wind  components.  Table  3.5a 
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Table  3.3  -  Results  from  optimum  weight  parameter  test  with 
DX  =  115  km. 


Root-Mean- Square-Error 


ght  Parameter 

Temperature 

U-  component 

V-  component 

12.5 

2.86  °F 

3.69  m/s 

3.12  m/s 

25. 

2.82 

3.55 

3.04 

50. 

2.79 

3.39 

2.90 

100. 

2.66 

3.16 

2.70 

250. 

2.89 

3.15 

2.50 

500. 

3.67 

3.47 

2.42 

1000. 

4.97 

3.96 

2.41 

1500. 

5.86 

4.33 

2.44 

2000. 

6.49 

4.68 

2.47 

3000. 

7.27 

4.97 

2.52 

5000. 

8.04 

5.34 

2.57 

Table  3.4  -  Same  as  Table  3.3  except  random  errors  have  been 
incorporated  into  the  observation  data,  with 
standard  deviation  of  random  error  of  2  '^F  and 
2  m/s  for  the  surface  temperature  and  wind 
components,  respectively. 


Root-Mean- Square- Err or 


ght  Parameter 

Temperature 

U-  component 

V-  component 

12.5 

3.05  °F 

3.61  m/s 

3.25  m/s 

25. 

2.82 

3.48 

3.06 

50. 

2.72 

3.32 

2.76 

100. 

2.51 

3.18 

2.45 

250. 

2.62 

3.23 

2.20 

500. 

3.40 

3.56 

2.19 

1000. 

4.75 

4.04 

2.26 

1500. 

5.70 

4.42 

2.32 

2000. 

6.36 

4.66 

2.37 

3000. 

7.19 

5.01 

2.44 

5000. 

8.09 

5.36 

2.57 
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much  less  than  the  u-  components. 

The  same  test  was  repeated  but  with  the  grid  spacing 

doubled.  The  purpose  was  to  determine  the  effect  on  the 

optimum  weight  parameter  and  the  least  RMSE  when  the 

interval  between  grid  points  is  changed.  The  results  in 

Table  3.3  show  that  the  optimum  k  for  each  variable 

o 

remains  the  same.  Also  worth  noting,  the  minimum  error  is 
only  slightly  less  when  using  the  denser  grid;  at  the  same 
time,  the  computation  required  is  four  times  greater,  which 
probably  does  not  justify  using  a  grid  spacing  much  less 
than  the  station  spacing,  in  this  case,  approximately 
95-100  km. 

3.3.2  Tests  on  data  fields  with  random  error 

To  substantiate  the  above  findings  on  optimum  weight 
parameters,  a  similar  test  is  conducted  with  additional 
random  error  added  to  the  fields.  The  random  errors  are 
generated  by  a  computer  and  have  a  Gaussian  distribution 
with  a  mean  of  zero  and  standard  deviation  of  2  °F  for 
temperature  and  2  m/s  for  the  wind  components.  This 
procedure  should  reveal  the  weight  parameter  which  yields 
the  optimum  smoothing  of  random  error,  that  is,  the  most 
effective  high  frequency  noise  filter.  The  actual  station 
values  without  the  random  errors  are  again  used  as  the 


standard . 
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Figure  3.4  - 


Graph  of  root-tnean-square-error  (RMSE)  for 
various  weight  parameters  (k  )  with  DX  =  57  km 
RMSEs  are  in  units  of  F  for  temp>erature  and 
m/s  for  wind  components.  (NOTE:  on  the  RMSE 
axis  the  linear  scale  changes  at  4.0) 


Table  3.2  -  Results  from  optimum  weight  parameter  test  with 
DX  =  57  km.  Root-mean-square-error  for  various 
weight  parameters  are  shown. 


Root-Mean- Square-Error 


Weight  Parameter  Temperature 


U-  component 


V-  component 
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I  for  this  station. 

3. )  Steps  1.)  and  2.)  are  repeated  for  each  of  a 
succession  of  missing  stations  chosen  at  random  from  the 

I  total  station  population. 

4. )  Then  an  RMSE  is  calculated  by  computing  the 
difference  of  the  observed  and  the  interpolated  values  at 

,  the  missing  station  for  each  analysis. 

5. )  Steps  1.)  through  4.)  are  then  repeated  for  each  of 

the  remaining  k  ' s  to  be  tested, 
o 

As  shown  in  Table  3.2  an  optimum  weight  with  a  minimum 
error  is  obtained  for  each  variable.  Fig.  3.4  shows  the 
graph  of  these  results;  rather  well-defined  minima  of  error 
j  are  depicted  in  the  case  of  temperature  and  u-  component  of 

the  wind.  For  temperature  a  of  approximately  100 
produced  the  minimum  error,  while  k^  lay  between  100  and 
j  250  for  u-  and  bet<«?een  500  and  600  for  the  v-  component  of 

the  wind.  These  results  suggest  that  a  weight  parameter  of 
approximately  250  would  be  optimal  for  the  meteorological 
variables  considered  in  this  particular  case.  This  type  of 
compromise  would  result  in  a  slight  "over-smoothing"  of  the 
temperature  field  and  "under-smoothing"  (noise)  of  the  v- 
component  wind  pattern.  The  apparent  discrepancy  of  the 
smaller  v-  field  RMSEs  as  com;-  '  d  to  those  of  the  u-  field 
is,  in  part,  due  to  the  fact  that  the  v-  components 
(generally,  the  east-west  components  since  the  y-  axis  of 
the  coordinate  system  is  defined  positive  to  the  east)  were 
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Barnes  analysis  is  performed  over  the  field  except  that  a 
single  station  observation  left  out.  Then  bilinear 
interpolation  of  the  four  surrounding  grid  values  from  the 
Barnes  analysis  is  used  to  obtain  an  estimate  at  the  missing 
station  location.  RMSEs  are  calculated  using  the  station's 
actual  observed  value  as  the  standard.  This  method  gives  a 
realistic  measure  of  error  since  the  datum  at  the  point  of 
interpolation  is  not  in  the  data  set  utilized  by  the 
analysis  scheme.  However,  if  this  datum  was  incorporated 
and  then  an  interpolation  made  back  to  the  point,  the  error 
could  approach  zero  as  approaches  zero  since  locally 
that  single  datum  almost  exclusively  influences  the 
surrounding  grid  point  values.  On  the  other  hand,  one  might 
expect  that  if  k^  becomes  too  small  in  a  "missing 
station"  test  short-wavelength  noise  will  increase  the 
error;  and  if  k^  becomes  too  large  excessive  smoothing 
will  occur  and  the  error  again  increases.  Therefore,  one 
can  reasonably  expect  an  optimum  k^  with  some  minimum 
RMSE. 

The  test  is  conducted  in  the  following  manner: 

1. )  Using  a  15X15  grid  (grid  spacing,  DX  of  57  km; 
one-half  of  that  used  in  the  "edge-effect"  test) 
encompassing  the  inner  35  stations  and  a  given  k^,  the 
objective  analysis  is  performed  with  one  station  missing. 

2. )  Through  bilinear  interpolation  of  the  four 
surrounding  grid  values,  an  estimated  value  is  calculated 
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response  of  the  features  in  the  analysis,  that  is,  the 
percentage  of  the  amplitude  of  each  wavelength  that  is 
retained.  The  weight  parameter  is  dependent  on  wavelength 
and  can  be  specified  for  the  smallest  resolvable  wave  (twice 
the  station  spacing  =  2  An)  when  the  response  for  the 
resolvable  scale  is  prescribed.  Barnes  (1973)  did  not  offer 
clear  guidance  for  the  selection  of  in  the  spatial 
weight  function.  His  choice  for  this  parameter  was  rather 
arbitrary  by  suggesting  that  one  should  know  beforehand  what 
the  response  of  the  smallest  resolvable  scale  should  be  as  a 
function  of  the  accuracy  or  representativeness  of  the  data, 

Koch  et  al.  (1981)  suggested  guidelines  for  choosing  the 
weight  parameter.  It  was  suggested  that  the  parameter  be 
selected  such  that  a  maximum  response  of  1/e  (.3679)  at  the 
2  An  scale  is  achieved.  However,  in  the  present  author’s 
opinion,  the  choice  of  1/e  as  the  maximum  response  for  the 
smallest  resolvable  wave  is  also  arbitrary.  Koch  did  not 
present  the  reasoning  behind  his  choice.  Therefore,  several 
tests  were  made  to  determine  an  optimum  weight  parameter  for 
a  given  data  distribution;  that  is,  one  that  yields  the 
least  error  in  the  analyzed  field. 

3.3.1  Optimum  weight  parameter  tests 

The  approach  taken  to  determine  an  optimum  weight 

parameter  (k  )  is  the  "missing  station"  method.  The 
o 


nonlinear  quasi-geostrophic  omega  and  vorticity  equations. 
Koelher  (1979)  also  used  Sanders'  model  with  certain 
modifications  in  order  to  provide  a  better  large-scale 
representation  of  meteorological  fields.  A  number  of 
modifications  are  employed  here  to  enhance  temperature 
gradients  and  to  sharpen  and  re-orient  height  (pressure) 
troughs, 

Sanders  defined  the  1000  mb  height  field  as  a  function 
of  X  and  y  (horizontal  distance  parameters) ,  and  the 
temperature  as  a  function  of  x,  y,  and  pressure.  His 
expressions  for  these  fields  are  given  in  Equations  4.1  and 
4.2,  where  x  and  y  in  his  original  functions  have  been 
replaced  by  their  latitudinal  and  longitudinal  equivalents. 

Z(X,4),1000)  »  2  cos  (  (X  -X+AX)  2it/L,  ]  COS  t  ( (^- 4>.^)  2ii/L  J  ]  (4.1) 

O  A  O  ^ 

T(X,4>,p)  =  jn<P)  “  tl-6ln  (1000/p)  ]  { nr  ( (^- (f.^)  + 

T  cos(  (X  -X)2ii/L,)  cost  («^-(t)„)2TT/L  .)  }  (4.2) 

ox  O  q) 

where  X  =  longitude  in  radians 
4)  =  latitude  in  radians 
p  =  pressure  in  millibars 

L  =  east-west  wavelength  in  radians  of  longitude 

A 

L  =  north-south  wavelength  in  radians  of  latitude 

4 

X^  =  reference  longitude 
4  =  reference  latitude 

AX  =  a  parameter  which  controls  the  phase  lag  of  the 
1000  mb  height  field  relative  to  the  temperature 
field  in  radians  of  longitude 


T^(p)  =  mean  temperature  at  pressure  p 
Z  =  amplitude  of  the  1000  mb  height  disturbance 
f  =  amplitude  of  the  1000  mb  temperature  field 
6  =  a  parameter  that  determines  the  pressure  at  which 
the  reversal  of  temperature  gradient  appears 
r  =  radius  of  the  earth 

ri  =  magnitude  of  the  mean  north-south  temperature 
gradient. 

Since  the  purpose  here  is  to  develop  a  two-dimensional 
model  of  meteorological  fields  at  or  near  the  surface  of  the 
earth,  the  temperature  expression  can  be  simplified  by 
letting  p  =  1000,  which  gives: 

T(X,^,1000)  *  T^(1000)  -  {nr((j)-(^^)  + 

T  cos  t  (X^-X)  )  cos  [  2tt/L^]  )  (4.3) 

This  function  proved  unsatisfactory  in  representing 
intense  thermal  gradients  in  frontal  zones;  a  nearly 
constant  gradient  over  the  field  results  when  using  Sanders' 
original  function.  With  coastal  fronts,  temperature 
gradients  of  5-10  °C  over  a  20  km  distance  are  not 
uncommon.  To  remedy  this  problem,  a  hyperbolic  tangent 
function,  dependent  on  a  temperature  difference,  is  utilized 
as  a  temperature  correction  which  is  added  to  the  pattern. 
The  correction  term  is  given  by: 

Tcop  =  B{TANH{  (T(X  ,((.,1000)-T^(1000)]/A))  (4.4) 

where  A  =  a  parameter  which  controls  the  intensity  of  the 
gradient  added  to  the  field 
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B  =  amplituc3e  of  the  temperature  correction. 

The  correction  function  is  designed  to  increase  temperatures 

greater  than  T  (1000)  and  decrease  temperatures  less 
m 

than  Tj^(1000).  Due  to  the  nature  of  the  hyperbolic 

tangent  function  which  has  its  first  derivative  maximizing 

at  zero,  the  temperature  gradient  (the  first  derivative  of 

temperature  with  respect  to  horizontal  distance)  becomes 

greatest  in  the  vicinity  of  the  chosen  isotherm  T  (1000) 

m 

where  T(^, 41, 1000)  -  T_(1000)  equals  zero. 

The  combined  expression  for  the  temperature  field  then 
becomes : 

=  T(X,4,1000)  +  T^qj^  .  (4.5) 

Sanders'  model  allowed  for  only  a  north-south  alignment 
of  troughs  and  ridges.  Therefore,  his  height  expression  had 
to  be  modified  to  obtain  any  desired  orientation  of  the 
trough  axis,  for  example,  along  the  coast  in  the  case  of  the 
coastal  front.  A  tilting  term  (Equation  4.6}  was 
incorporated  which  forces  the  longitude  to  vary  as  a 
function  of  latitude  difference  from  the  reference  latitude: 


X  =  c  ( 4—  4  ) 

TILT  '  o' 


(4.6) 


where  c  is  a  parameter  which  controls  the  degree  of  tilting, 
that  is,  the  slope.  Sanders'  phase  lag  parameter  (^^)  is 
not  used  here  since  separate  reference  latitudes  and 


longitudes  are  defined  for  the  temperature  and  height 


fields. 

To  make  the  height  gradient  pattern  around  maxima  and 


minima  (highs  and  lows)  more  realistic  (by  increasing  the 

gradient  north  and  decreasing  the  gradient  south  of  the 

height  center),  a  linear  correction  is  added  to  the 

latitudinal  factor  of  the  height  expression,  where 

cos  (  (  )  2ti/L  .  ] 

o  c 

is  replaced  with 

cos  ( (<{>-<f^)  2Tr/L^]  +  D('(>-4'  ) 

and  D  is  a  parameter  which  controls  the  degree  of  change 
applied  to  the  height  gradient. 

Finally,  a  height  correction  dependent  on  the  absolute 
difference  of  temperature  from  a  given  point  to  a  chosen 
isotherm  is  added  to  the  height  field.  The  chosen  isotherm 
is  preferably  one  located  in  the  frontal  zone.  Thus,  the 
smallest  corrections  are  being  added  near  the  front.  This, 
in  effect,  forces  the  trough  to  be  sharpened  in  the 
immediate  vicinity  of  the  front  assuming  that  the  patterns 
are  matched  beforehand  with  the  temperature  gradient  lying 
in  the  trough  of  the  height  field.  The  height  correction 
due  to  temperature  is  a  quadratic  function  given  by: 

HTCOR  =  (4.7 

where  E  =  a  parameter  which  controls  the  degree  of 
sharpening  of  the  trough,  and 
TpRObjT  ®  mean  temperature  (isotherm)  chosen  to 

represent  the  frontal  zone. 

The  final  expression  for  the  height  field  thus  becomes: 


Z  ,1000  ) 


=  Z  COs{  +  X 

{cost  )2n/L.]+D(C-<  ) }  +  HTCOR.  (4.8) 

Sample  plots  of  the  resultant  temperature  and  height  fields 
are  shown  in  Fig.  4.1.  One  should  note  that  the  surface 
trough  through  Missouri  is  at  the  warm  edge  of  the 
temperature  gradient;  that  is,  along  the  surface  cold  front. 
Also  note  the  production  of  a  warm  frontal  trough  through 
Michigan  and  Ohio. 

Due  to  the  modifications  mentioned  above  the  functions 
of  temperature  and  height  have  become  increasingly  complex. 
Temperature  is  now  a  composite  function  with  the  hyperbolic 
tangent  incorporated,  and  the  height  field  is  dependent  on 
this  temperature.  Therefore,  finite-differencing  will  be 
used  to  create  the  simulated  fields  of  temperature  gradients 
and  of  the  wind.  Finite-differencing  of  the  resulting  wind 
components  will  be  used  to  obtain  additional  kinematic 
variables  such  as  divergence,  vorticity,  and  temperature 
advection . 


4.2  The  Simulated  Wind  Field 


The  geostrophic  wind  ,  which  is  derived  directly 
from  the  height  field,  is  the  basis  for  the  simulated  wind 
field  V--...  The  geostrophic  wind  components  are  defined 

Sin 

by: 


V 

g 


(-g/f)  (''Z/.‘y)  ; 


(g/f) (oz/Lx) 


(4.9) 
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Figure 


Example.of  simulated  1000  mb  temperature 
f ield  (  F)  . 


Figure  4.1b  -  Example  of  simulated  1000  mb  height  field  (m)  . 
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where  g  is  the  acceleration  of  gravity 
f  is  the  Coriolis  parameter 

z  is  the  1000  mb  height  from  the  analytic  field. 
Friction  is  then  added  fictitiously  to  the  geostrophic 
wind  by  utilizing  a  linear  relation  between  the  "actual" 
wind  speed  and  the  geostrophic  wind  speed,  and  by  reducing 
the  geostrophic  wind  direction  by  a  constant  amount  (a)  and, 
thus,  creating  cross-isobar ic  flow  toward  low  pressure. 

These  alterations  on  the  geostrophic  wind  are  as  follows: 


*^SIM*DIR  ~  *^g^DIR  ‘  “ 


if  •''g’DIR  -  “  >  ^ 


'  <%>DIR  -  “>♦36®°  if  (^g'DIR  -  «  <  ».  (4.11) 

Baur  and  Phillips  (1938)  computed  the  angle  (a)  between 
the  surface  wind  and  the  geostrophic  wind  (land  and  sea) , 
and  the  ratio  (a)  of  the  surface  wind  speed  and  the 
geostrophic  wind  speed  for  different  latitudes  (Table  4.1). 
For  the  test  cases  in  this  research,  mean  values  (those  for 
the  central  latitude  of  the  grid)  were  chosen  for  all  wind 
computations.  For  example,  in  the  case  of  the  simulated 
coastal  front  in  the  PAM  II  network,  values  of  43,  17,  .35, 
and  .60  were  used  for  a  ,  a  ,  a_ ,  and  a  ,  respectively. 

A  review  of  papers  on  east  coast  f rontogenesis  by  Bosart  et 
al.  (1972)  and  Bosart  (1975)  confirmed  these  rough 
approximations. 


The  other  contribution  to  the  ageostrophic  component  of 
the  wind  which  results  from  acceleration  was  also  considered 
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Table  4.1  -  Th^  angle  (a)  between  the  surface  wind 
wind  and  the  geostrophic  wind,  and  the 
ratio  (a)  between  the  surface  wind 
speed  and  the  geostrophic  wind  speed 
for  different  latitudes  ('J'),  according 
to  Baur  and  Philipps  (1938,  p.  292). 


“l 

“s 

a 

s 

0 

61 

40 

•  •  •  • 

•  •  •  • 

10 

55 

29 

•  •  •  • 

•  •  •  • 

20 

49.5 

23 

•  •  •  • 

•  •  •  • 

30 

45 

19 

0.31 

0.56 

40 

42 

16 

0.38 

0.63 

50 

39 

14.5 

0.42 

0.67 

60 

37 

13 

0.46 

0.70 

70 

36 

12.5 

0.485 

0.715 

80 

35 

12 

0.495 

0.723 

90 

35 

12 
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for  inclusion  in  the  generation  of  a  realistic  surface  wind. 

The  ageostrophic  wind  V  due  to  this  cause  was  derived 

d 

as  follows: 

dv/dt  =  V  ‘VV  (4.12) 

g  g 

assuming  the  local  time  rate  change  of  momentum  and  the 
vertical  advection  of  momentum  are  zero,  and  replacing  V 
with  in  the  horizontal  advection  term.  The 

g 

acceleration  is  also  given  by: 

dVdt  =  -flc  X  (V-V  )  =  -fie  XV  ;  (4.13) 

g  a 

therefore , 

v.^  =  H/f  X  (Vg-^Vg)  .  (4.14) 

An  experiment  was  conducted  to  compare  the  effects  of 
the  fictitious  "friction"  and  the  ageostrophic  wind  0 

d 

on  the  divergence  field.  The  geostrophic  wind  field  was 
derived  from  the  1000  mb  height  pattern  shown  in  Fig.  4.1b. 
Then  Baur  and  Phillips'  frictional  modification  and  V 

d 

were  added  separately  to  this  basis  geostrophic  wind  field. 
From  these  resulting  fields  divergences  were  obtained.  A 
comparison  of  grid  point  values  revealed  that  the 
contribution  to  the  divergence  field  by  V  was  small 
compared  to  the  contribution  due  to  "friction". 

Additional  support  for  not  incorporating  d^/dt  comes 
from  a  comparison  of  the  relative  magnitudes  of  the  forces 
acting  on  the  wind.  The  Coriolis  force  ,CF=f|V|,  based  upon 
the  frictionless  wind,  V-Vg+V^;  the  horizontal  pressure 
gradient  force,  HPGF=-gvz;  and  the  accelerational  force  due 
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to  the  advection  of  the  geostrophic  wind  (Equation  4.12) 
were  computed  over  a  15X18  grid.  The  Coriolis  force  and  the 
horizontal  pressure  gradient  force  were  an  order  of 
magnitude  greater  than  d\^/dt  (from  Equation  4.12)  over 
approximately  80-90%  of  the  field.  The  accelerational  force 
exceeded  each  of  the  other  two  forces  at  only  2  of  the  154 
test  points.  Furthermore,  the  wind  field  resulting  from  the 
approximation  given  by  Baur  and  Phillips  exhibits  a  very 
reasonable  pattern  and  is  more  than  sufficient  for  this 
study. 

4.3  Sensitivity  of  Analyzed  Fields  to  Error  Sources 

The  simulated  fields  of  temperature  and  height  described 
above  were  used  to  compare  the  error  in  analyzed  fields  that 
would  result  from  random  error  and  from  spatial  data 
distribution  and  interpolation  by  the  objective  analysis 
scheme.  The  analytic  meteorological  fields  provide  a  known 
standard  for  each  point  of  the  grid.  To  accomplish  the 
sensitivity  test  the  following  procedure  was  employed: 

1. )  The  Barnes  analysis  was  performed  on  the  bogus 
station  data  without  random  error. 

2. )  The  analysis  was  repeated  with  random  error 
superimposed  on  the  data  set.  Random  errors  with  mean  of 
zero  and  standard  deviation  of  2  °F  and  3  m  were  added 

to  all  of  the  simulated  observation  values  of  temperature 


and  height,  respectively. 
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3.)  The  two  analyses  were  compared  to  the  "actual" 
simulated  grid  point  values,  the  standard,  in  order  to 
calculate  RMSEs. 

Table  4.2  shows  that  the  error  in  temperature  increased 
slightly,  while  the  error  in  height  actually  decreased  when 
random  errors  were  introduced  into  the  fields.  As  the 
latter  result  was  unexpected  a  second  test  was  conducted 
with  a  larger  standard  error  of  height,  5  m.  This  led  to  a 
very  slight  increase  in  mean  error.  These  results  reveal 
that  the  analysis  scheme  smoothes  random  error  very 
efficiently.  That  is,  the  error  due  to  random  error  is 
practically  nil  compared  to  the  error  resulting  from  the 
data  distribution  and  the  analysis  scheme's  interpolation 


procedure 


Table  4.2  -  Analysis  error  due  to  random  error.  Root-mean- 
square-errors  are  shown  for  analyses  with  and 
without  random  error  using  simulated  fields  of 
surface  temperature  and  1000  mb  height. 


Temperature  1000  mb  Height 


Without  Random  Error 


1.37  °F 


4.74  m 


With  Random  Error-- 

(Standard  deviation  of  1.73  4.66 

error:  2  °F  and  3  m) 


(Standard  deviation  of 
error:  5  m) 


4.75 
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CHAPTER  5 

SENSITIVITY  OF  ANALYZED  METEOROLOGICAL  FIELDS 
TO  OBSERVATIONAL  NETWORK  DESIGN 

5.1  Simulated  "True"  Meteorological  Fields 

The  mathematical  functions  defined  in  the  previous 
chapter  were  used  to  produce  meteorological  fields 
representative  of  a  coastal  front  typical  of  the  United 
States  east  coast  in  winter.  A  theoretical  measure  of  the 
analysis  accuracy  of  different  observing  systems  can  be  made 
since  the  "true"  values  of  all  meteorological  variables  are 
known  over  the  entire  field.  The  purpose  here  is  to 
determine  how  well  the  proposed  PAM  II  network  for  the 
winter  1986  GALE  project,  in  combination  with  the  Barnes 
objective  analysis  scheme,  describes  this  mesoscale  weather 
system.  Comparisons  are  made  with  four  alternative 
observing  networks  to  gain  insight  into  possible 
improvements  or  changes  to  the  proposed  network  design. 
Coastal  fronts  are  investigated  since  a  large  portion  of 
GALE  will  concern  the  study  of  these  mesoscale  systems  and 
their  role  in  cyclogenesis.  Coastal  fronts  provide  intense 


5E 

temperature  gradients  and  strong  horizontal  wind  shears 
which  offer  a  means  to  fully  test  the  observing  network's 
ability  in  describing  related  "true"  meteorological  fields. 

Fig.  5.1  depicts  the  simulation  of  an  onshore  front  over 
the  coastal  plains  of  North  Carolina  and  South  Carolina.  A 
realistic  coastal  front  situation  is  duplicated  with  a  wedge 
of  cold  air  between  the  front  and  the  Appalachian  Mountains; 
the  cold  air  wedge  is  associated  with  the  ridge  of  high 
pressure.  The  coastal  front  lies  in  a  surface  trough  and  is 
squeezed  between  this  cold  air  dome  and  the  warmer  air  over 
the  ocean  resulting  in  a  temperature  gradient  in  excess  of 
4  °C  over  10  km. 

These  patterns  were  produced  using  values  (obtained  from 
the  analytic  functions)  on  a  regular  rectangular  grid  point 
array  with  grid  spacing  of  12.5  km.  The  grid  is  positioned 
on  a  coordinate  system  which  has  its  x-axis  perpendicular, 
and  y-axis  parallel,  to  the  coast.  This  allows  the  grid 
array  to  be  aligned  with  the  proposed  PAM  II  network  which 
is  described  in  the  next  section. 

The  "true"  kinematic  fields,  derived  from  the  simulated 
wind  field  (Fig.  5.1c),  and  the  temperature  gradient  field 
(magnitude)  are  depicted  in  Fig.  5.2.  These  patterns  are 
representive  of  an  intense  coastal  front.  An  axis  of 
maximum  positive  vorticity,  with  values  in  excess  of 
30  X  10”^/s,  is  aligned  along  the  frontal  zone  over  the 
Carolines.  Equally  strong  convergence  is  depicted  over  this 
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Figure  5.6a  -  Barnes  analysis  of  temperature  field  (°C)  for 
simulated  coastal  front  case. 


'igure  5.6b  -  Barnes  analysis  of  1000  mb  height  field  (rr) 
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Table  5.2  -  Effects  of  random  error.  (Same  as  Table  5.1 
except  random  error  is  added  to  PAM  Il^data, 
with  standard  deviation  of  error  of  2  C,  8  m, 
and  2  m/s  for  temperature,  1000  mb  height,  and 
wind  components,  respectively.) 


RMSE 

P.  E 

Temperature 

2.57 

1000  mb  Height 

3.91 

U-  component 

1.97 

36.1 

V-  component 

1.85 

37.3 

Divergence 

4,53 

76.9 

Pel .  Vort icity 

5.73 

66.7 

Temp.  Gradient 

0.67 

80.6 

Temp.  Advection 

3.63 

109.8 

Local  Frontogenes i s 

1.11 

142.5 
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analyses  confirmed  the  previous  results  in  Chapters  3  and  4. 
Random  errors  had  little  effect  on  the  RMSE  as  shown  in 
Table  5.2;  that  is,  the  Barnes  scheme  effectively  filtered 
out  these  errors  in  the  data.  For  example,  the  RMSE  for  the 
temperature  field  increased  only  0.17  °C,  from  2.40  °C 
(without  random  error)  to  2.57  °C  (with  random  error).  The 
P.  E.’s  increased  only  1-6%  in  the  kinematic  parameters. 

The  results  in  Table  5.2  reveal  the  magnitude  of  typical 
errors  to  be  expected  with  a  network  such  as  PAM  II.  If  the 
RMSE  is  compared  to  the  average  magnitude  of  the  parameter 
over  the  whole  objectively  analyzed  field,  the  percentage 
error  would  be  unacceptably  large,  100-300%.  That  is,  the 
error  would  be  1  to  3  times  as  large  as  the  analyzed  value. 

Fig.  5.6  shows  analyses  resulting  from  using  the 
proposed  PAM  II  network  data  with  random  errors  incor¬ 
porated.  These  analyses  should  be  compared  to  Figs.  5.1 
and  5.4.  It  is  obvious  that  the  network  does  not  adequately 
represent  the  true  meteorological  fields.  The  Barnes 
analyzed  temperature  and  1000  mb  height  fields  are  grossly 
over-smoothed  (compare  Figs.  5.1a,b  with  Figs.  5.6a,b).  The 
same  is  true  in  the  cases  of  relative  vorticity,  divergence, 
and  temperature  gradient  (Figs.  5.6d-f),  for  which  the 
respective  maxima  are  6  X  10~^/s,  -4  X  10  ^/s,  and 
0.75  °C/10  km.  These  values  are  only  about  30%  of  the 
magnitudes  of  these  maxima  in  the  true  fields. 
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Table  5.1  -  Barnes  analysis  (using  simulated  PAM  II  data)  vs. 

"true"  field.  Root-mean-square-error  (RMSE)  and 
percentage  error  (P.  E.)  are  shown. 


RMSE 

P.  E 

Temperature  (°C) 

2.40 

1000  mb  Height  (m) 

2.71 

U-  component  (m/s) 

1.93 

35.3 

V-  component  (m/s) 

1.67 

33.6 

Divergence  (10”^/s) 

4.28 

72.6 

Rel .  Vorticity  (10~^/s) 

5.22 

60.8 

Temp.  Gradient  (°C/10^m) 

0.64 

77.7 

o  4 

Temp.  Advection  (  C/10  s) 

3.57 

108.1 

o  8 

Local  Frontogenes i s  (  C/10  ms) 

1.10 

140.5 
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described  above.  However,  since  the  "true"  field,  which 
serves  as  the  standard,  is  represented  by  values  on  the 
finite  grid,  the  simulated  station  values  in  the  present 
experiments  are  obtained  from  this  same  set  of  gridded 
values  using  bilinear  interpolation. 

These  PAM  II  station  "data"  are  now  objectively 
analyzed.  Table  5.1  summarizes  the  RMSEs  and  percentage 
errors  of  the  analysis  based  on  the  proposed  PAM  II  network 
compared  to  the  true  field  at  grid  points.  The  percentage 
error  is  given  by: 

M 

P.  E.  =  RMSE/( (  L  IV^P/M)  (5.1) 

m=l 

and  V  is  the  grid  point  value  from  the  true 
m 

meteorological  field  and  M  is  the  number  of  grid  points. 

The  RMSE  of  divergence  over  the  proposed  PAM  II  region  is 
about  73%  of  the  average  magnitude  of  the  divergences  over 
the  true  field  and  similarily,  about  61%  for  the  relative 
vorticity.  For  temperature  advection  and  local 
f rontogenesis  (local  time  rate  of  change  of  the  magnitude  of 
the  temperature  gradient)  the  P.  E.  exceeds  100%. 

Random  errors  with  means  of  zero  and  standard  deviation 
of  error  of  2  °C,  8  m,  and  2  m/s  for  temperature, 
height,  and  u-,  v-  wind  components,  respectively,  were 
incorporated  into  the  PAM  II  station  values.  Such  errors 
might  simulate  typical  instrumental  and  calibration  errors, 
as  well  as  local  exposure  and  siting  effects.  The  subsequent 
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Figure  5.5  -  Proposed  42-’Station  PAM  II  network  with 

approximately  60  km  spacing  normal  to  the 
coast  and  75  km  parallel  to  the  coast. 


T 
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dense  array  from  the  true  field  (DX  =  12.5  km).  Therefore, 
in  order  to  make  comparisons  between  the  "true"  fields 
(which  serve  as  the  standard)  and  the  objectively  analyzed 
data  as  fair  as  possible,  the  temperature  and  wind 
components  from  the  true  fields  are  f i n i te-d i f ferenced  using 
a  DX  of  25  km,  that  is,  using  every  other  grid  point.  This 
naturally  decreases  the  magnitude  of  the  extreme  values  of 
the  various  parameters  in  the  vicinity  of  the  front  and 
along  the  coast.  The  patterns  shown  in  Fig.  5.4  are  used  as 
the  standard  against  which  objectively  analyzed  fields  are 
compared.  The  maximum  values  in  these  fields;  20  X  10~^/s 
for  vorticity,  -15  X  10  ^/s  for  divergence,  and  2.5  °C/10  km 
for  temperature  gradient;  are  only  50-60%  as  large  as  the 
maxima  in  the  fields  where  DX  *  12.5  km  (Fig.  5.2). 

5.2  Experiments  with  the  Proposed  PAM  II  Network 

The  proposed  PAM  II  network  of  GALE  shown  in  Fig.  5.5  is 
made  up  of  42  land-based  stations  in  a  near-rectangular 
array.  The  values  of  meteorological  variables  (temperature, 
1000  mb  height,  and  wind  components)  at  these  stations  are 
obtained  by  bilinearly  interpolating  the  four  surrounding 
grid  point  values  from  the  true  simulated  field  where 
DX  =  12.5  km.  The  station  values  were  also  obtained 
directly  using  the  mathematical  functions.  This  gave  values 


essentially  equal  to  those  obtained  using  the  procedure 
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Figure 


5.3  -  Example  of  cellular  pattern  which  results  from 
using  a  finite  number  of  grid  values.  Dashed 
lines  are  contours  drawn  when  considering  only 
grid  point  values.  Solid  lines  represent  the 
"true"  continuum  field. 
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region  in  Fig.  5.2b,  The  magnitude  of  the  temperature 
gradient  exceeds  4  °C/10  km  across  the  front  (Fig.  5.2c). 

The  cellular  patterns  in  these  fields  arise  from  the 
inherent  problem  of  trying  to  represent  a  continuous  field 
with  a  finite  array  of  grid  points.  The  alternative  maxima 
and  minima  in  the  pattern  produced  by  the  contouring  scheme 
(SURFACE  II)  are,  as  they  must  be,  particularly  prominent 
along  the  frontal  zone  where  gradients  are  strongest.  This 
problem  of  cellular  patterns  can  be  demonstrated  with  a 
simple  example  as  shown  in  Fig.  5.3. 

If  the  grid  network  was  dense  enough,  that  is,  DX 
approaches  zero,  the  contours  would  appear  as  the  solid 
lines  in  the  figure;  however,  since  only  a  finite  number  of 
grid  points  is  available  the  true  axis  of  maximum  values 
falls  between  the  grids  points  and  cannot  be  correctly 
represented.  The  result  is  a  cellular  pattern  of  maxima 
along  this  axis.  To  verify  that  this  feature  does  not  exist 
in  the  true  continuum,  the  fields  were  reconstructed 
allowing  the  front  to  parallel  the  grid  rows;  with  this 
orientation  the  cellular  pattern  was  removed.  In  the 
present  experiments  the  gridded  fields  with  the  cells  are 
accepted  as  the  "true"  fields  with  the  knowledge  that  the 
atmosphere  is  not  quite  correctly  represented. 

In  the  tests  to  be  described,  station  values  are 
interpolated  to  points  of  a  grid  which  has  a  grid  spacing  of 
25  km.  This  grid  array  is  simply  a  subset  set  of  the  more 


Figure 


Figure  5.1a  -  Simulated  temperature  field  (°C)  depicting  an 
onshore  coastal  front. 
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Figuie  5.6c 


Vector  display  of  Barnes  analysis  of  wind 
field. 


Figure  5.6d  - 


Vertical  component  of  relative  vorticity 
(10  /s)  derived  from  Barnes  analyzed  wind 

component  s . 


Divergence  field  (10  /s)  derived  from  Barnes 

analyzed  wind  components. 


Temperature  gradient  field  {°C/10  km)  derived 
from  Barnes  analyzed  temperature  field. 


7  5 

5.3  Alternative  PAM  II  Networks 

Is  the  proposed  PAM  II  network  the  configuration  which 
gives  the  best  results  as  regards  reasonable  resolution  of 
coastal  front  features?  To  answer  this  question  three 
alternative  networks  covering  the  same  geographic  area  are 
investigated : 

1. )  9X5  network  -  is  a  network  configuration  (shown  in 
Fig.  5.7a)  which  is  simply  a  re-arrangement  of  the  42  PAM  II 
stations  enhancing  resolution  normal  to  the  coast  (spacing 
of  30  km)  at  the  expense  of  that  along-coast  (spacing  of 
150  km) ; 

2. )  9  X  17  network  -  has  approximately  four  times  the 
number  of  stations  as  in  the  present  PAM  II  network;  that 
is,  the  station  density  is  doubled  in  both  directions 
resulting  in  station  spacings  of  30  km  normal  to  the  coast 
and  37.5  km  parallel  to  the  coast  (Fig.  5.7b); 

3. )  5X5  network  -  contains  approximately  one-half  the 
number  of  PAM  II  stations  (spacing  of  60  km  normal  to  and 
150  km  parallel  to  the  coast)  and  is  more  nearly  compatible 
with  the  spacing  of  the  regular  hourly  reporting  stations 
(Fig .  5.7c). 

In  addition,  a  network  of  34  regular  hourly  reporting 
stations  (National  Weather  Service,  Federal  Aviation 
Administration,  and  military  stations)  overlying 
approximately  the  same  area  was  subjected  to  the  same  test 


Figure  5.7a  -  A  possible  ”9  X  5”  PAM  II  observation  network. 


Figure  5.7d  -  Regular  hourly  reporting  station  network. 
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(Fig.  5.7d)  . 

These  tests  were  made  with  random  error  added  to  the 
station  values.  A  range  of  weight  parameters  ,k^'s, 
were  used  in  the  tests  because  the  different  station 
densities  warranted  it.  However,  to  make  equitable 
comparisons  between  networks,  results  using  the  same  k^ 
are  also  compared. 

A  summary  of  the  results  of  these  several  experiments 

are  outlined  in  Table  5.3.  It  is  apparent  that  the  RMSE 

decreased  with  increasing  station  density.  The  least  dense 

network,  5X5  network  (23  stations)  produced  the  largest 

errors  with  a  mean  error  of  17.41  m  for  the  1000  mb  height 

and  12.60  m/s  for  the  u-  component  of  the  wind  (k^  =  200). 

The  9  X  17  configuration  generally  produced  the  smallest 

errors  which  were  about  10%  less  than  those  generated  using 

the  proposed  PAM  II  (5  X  9)  network.  For  example,  in  the 

relative  vorticity  field  an  error  of  5.11  X  10  ^/s 

-5 

occurred  with  the  9  X  17  network  and  5.73  X  10  /s  with 
the  proposed  PAM  II.  However,  this  is  only  a  slight 
improvement  when  one  considers  that  the  station  density  was 
approximately  quadrupled.  Since  the  scale  of  the  mesoscale 
features  of  the  front  are  small  compared  to  the  average 
station  spacing,  the  9  X  17  network  fails  to  significantly 
improve  the  analysis  accuracy. 

The  re-arrangement  of  the  proposed  PAM  II  network  into 
the  9X5  configuration  increased  the  station  density  across 
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the  front  and  decreased  the  density  along  the  frontal 
boundary.  As  a  result,  the  analysis  accuracy  changed 
little;  the  mean  errors  approximately  equaled  those  of  the 
PAM  II  network  (Table  5.3).  While  offering  a  potential 
monetary  savings  through  requiring  fewer  observational 
stations,  the  5X5  network  proved  to  be  significantly  less 
accurate  than  the  proposed  PAM  II  network. 

The  proposed  PAM  II  system  was  more  accurate  than  the 
regular  hourly  reporting  network.  This  is  partly  due  to  the 
greater  station  density  of  the  PAM  II  network  (42  vs.  34 
stations),  and  partly  to  the  irregular  spatial  distribution 
of  the  hourly  reporting  stations.  The  PAM  II  system  has 
another  important  advantage  over  the  regular  reporting 
stations  which  was  not  simulated  in  this  study.  PAM  II 
stations  will  be  able  to  make  observations  at  small  time 
intervals  (every  five  minutes) .  This  increased  time 
resolution  will  provide  data  which  can  be  filtered  with 
respect  to  time  to  remove  high  frequency  noise  and  make  the 
analysis  more  reliable.  The  regular  reporting  stations  only 
make  one  observation  per  hour  which,  in  most  cases,  will  not 
be  representative  of  the  whole  hour.  Furthermore,  the 
greater  time  resolution  of  the  PAM  II  network  will  permit 
t ime- to- space  conversion  of  observations  and  thus  enhance 
spatial  resolution  if  features  move  with  known  velocity 
through  the  observation  network. 

The  results  in  Table  5.3  show  that  the  mean  error 


6  3 

generally  decreased  as  the  weight  parameter  decreased.  This 
is  expected  when  error-free  data  taken  from  the  true  field 
are  analyzed.  However,  even  with  random  errors  in  the  data 
the  RMSE  can  continue  to  decrease  with  decreasing  k 

o 

because  the  component  of  the  RMSE  independent  of  random 

error  dominates  over  the  component  of  the  RMSE  due  to  random 

error.  The  former  component,  which  is  due  to  the  nearly 

unresolvable  short-wave  features  in  the  frontal  zone,  is  so 

large  that  the  incorporation  of  random  errors  has  little 

effect  on  the  noise  level  of  the  field.  Therefore, 

increased  noise  in  the  field  due  to  random  error  does  not 

become  the  dominate  feature  and  allow  an  optimum  k  to 

o 

be  obtained,  as  was  the  case  with  the  optimum  weight 
parameter  test  in  Chapter  3. 

The  above  experiments  were  repeated  for  each  station 
network  without  incorporating  random  error.  Table  5.4  shows 
RMSEs  calculated  when  these  analyses  are  compared  to  the 
true  fields.  The  analyses  of  the  uncontaminated  data  fields 
lead  to  results  similar  to  those  in  Table  5.3.  Again,  the 
mean  error  is  correlated  to  station  density,  with  the 
proposed  PAM  II  network  having  errors  less  than  the  5X5 
and  hourly  reporting  station  networks,  approximately  equal 
to  the  9X5  network,  and  slightly  greater  than  the  9X1'’ 
network . 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 

The  process  of  objective  analysis  has  been  investigated 
in  this  study  using  the  interpolation  technique  devised  by 
Barnes  (1973).  A  method  of  determining  optimum  weight 
parameters  for  the  Barnes  objective  analysis  scheme  has  been 
tested.  The  "missing  station"  technique  enabled  a  single 
value  of  weight  parameter  to  yield  a  minimum  mean  error  in 
the  analysis  since  the  datum  at  the  point  of  interpolation 
was  not  in  the  data  set  utilized  by  the  analysis  scheme. 

This  method  seems  more  appropriate  than  those  employed  by 
Barnes  (1973)  and  Koch  et  al.  (1981).  The  tests  described 
in  Chapter  3  showed  that  the  weight  parameter  that  resulted 
in  a  total  response  of  0.936  for  the  minimum  resolvable  2  tn 
wave  yielded  the  optimum  results.  However,  these  tests  were 
performed  with  only  one  case;  several  are  needed  to  obtain 
more  conclusive  results.  Therefore,  it  is  suggested  that  a 
test  similar  to  the  one  described  be  conducted  prior  to 
extensive  analyses  to  obtain  the  appropriate  weight 
parameter  for  a  given  data  distribution  and  meteorological 


setting. 


Analyses  produced  on  a  arid  mesh  with  grid  spacing 
approximately  equal  to  the  average  station  spacing  were 
compatible  with  those  produced  on  the  grid  with  spacing  of 
about  one-half  the  initial  grid.  Therefore,  when  costs  are 
considered,  it  is  more  practical  to  use  a  DX  >  Ln/2. 

The  Barnes  objective  analysis  scheme  proved  generally 
more  accurate  than  detailed  hand  analyses  of  kinematic 
fields  such  as  divergence  and  vorticity.  These  comparisons 
were  made  using  the  average  of  several  subjective  analyses 
performed  by  experienced  analysts  as  the  standard. 

Edge  effects  are  quite  pronounced;  that  is,  the  Barnes 
interpolation  became  unreliable  near  the  edge  of  the  data 
domain.  Extrapolation  into  data  void  regions  gave 
increasingly  erroneous  results.  This  mimics  similar 
difficulties  in  subjective  analyses. 

When  data  sets  were  artificially  contaminated  with 
random  errors  of  a  magnitude  which  simulates  typical 
instrument  and  exposure  errors,  the  final  analysis  revealed 
no  deterimental  effects.  That  is,  the  Barnes  scheme 
effectively  smoothed  (filtered)  the  random  error  from  the 
analysis.  As  a  result,  most  of  the  error  in  the  analysis 
can  be  attributed  to  the  application  of  the  analysis  scheme 
(interpolation)  and  to  the  limitations  of  the  observing 
network,  for  example,  its  configuration,  station  density, 
and  frequency  of  observation.  However,  different  results 
could  be  obtained  if  the  random  errors  added  to  the  field 
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are  much  larger  than  values  used  in  this  study.  After  the 
incorporation  of  an  extremely  high  degree  of  artificial 
error,  the  process  of  removing  that  error  would  inevitably 
remove  much  of  the  texture  of  the  original  data  field. 

A  mathematical  model  has  been  developed  which  enables 
fronts  with  extremely  intense  gradients  to  be  constructed. 
Through  the  simulation  of  a  typical  east  coast  coastal 
front,  the  analysis  accuracies  of  the  proposed  (5  X  9) 

PAM  II  network  for  the  1986  GALE  project  and  four 
alternative  networks  were  tested.  In  addition,  an  a  priori 
estimate  of  typical  errors,  which  results  from  utilizing  the 
PAM  II  network  and  the  Barnes  objective  analysis  as  a  tool 
to  reconstruct  meteorological  fields,  has  been  accomplished. 

These  experiments  show  that  errors  are  unacceptably 
large  in  the  kinematic  fields  on  the  mesoscale,  with  RMSEs 
approaching  100%  of  the  true  values,  and  even  greater  in 
temperature  advection  and  local  frontogenesis.  Therefore, 
one  is  fortunate  just  to  achieve  the  correct  sign  (+  or  -) 
of  the  divergence  and  vorticity  when  attempting  to 
reconstruct  the  mesoscale  pattern  with  the  present  PAM  II 
design . 

Though  being  poor  in  its  description  of  the  "true” 
fields,  the  proposed  PAM  II  network  is  nearly  as  good  as  a 
9  X  17  configuration  which  halves  the  station  spacing.  The 
more  dense  9  X  17  network  did  not  provide  a  significant 
improvement,  probably  because  the  scale  of  the  features  in 
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the  fields  was  so  small  that  the  increase  in  station  density 
had  little  effect;  that  is,  the  stations  were  still  too  far 
apart  to  adequately  resolve  the  intense  gradients.  The 
PAM  II  network  (42  stations)  produced  somewhat  better 
results  than  the  regular  reporting  hourly  network  (34 
stations).  However,  the  potential  exists  that  the  analysis 
accuracy  of  the  proposed  PAM  II  network  can  be  significantly 
better  than  that  of  the  regular  station  network  since  the 
PAM  II  stations  will  be  able  to  make,  record,  and  transmit 
observations  almost  continuously.  Time-to-space  conversion 
of  observations  can  be  helpful  if  the  coastal  front,  or  any 
feature,  is  moving  and  is  approximately  in  steady-state. 

This  would,  in  effect,  increase  the  data  density 
considerably.  Also  automobile  transverses  across  the  front, 
not  simulated  in  this  study,  could  possibly  improve  the 
analysis  accuracy. 

In  order  to  fully  depict  mesoscale  features,  the  results 
of  experiments  in  this  research  indicate  that  spatial 
separation  of  observation  sites  should  be  on  the  order  of 
the  wavelength  (or  less)  of  the  mesoscale  feature  under 
consideration.  That  is,  for  an  observational  network  to 
describe  the  intense  temperature  gradients  and  wind  shears 
typical  of  a  coastal  front  and  which  occur  in  a  narrow  zone 
of  5  to  20  km  in  width,  observations  are  needed 
approximately  every  5  to  10  km,  or  even  less. 
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